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Capitolul | 
_ ERORI. CONCEPTE DE BAZĂ IN CALCULUL NUMERIC 


_$ 1. Introducere 


_ Calculul aproximativ joacă un rol central în matematica aplicată. 
_ În practică ne mulțumim foarte adesea cu o valoare aproximativă a can- 
ăţii pe „care o găsim“, alteori ienorăm unele cantități mici în comparaţie cu 
tele, iar nu de puţine ori „datele de intrare“, ale problemei puse, sînt inexacte. 
- Introducem următoarele notații (nu sînt exacte din punct de vedere 
atematic) utile în practică: 

a < b (sau d > a) citim: „a este mult mai mic decit b* (sau „b este mult 
i mare decit a”). - 
Evident, acest „mult mai mare“ sau „mult mai mic“ depinde de contextul 


problemei (mneori ab dacă a< - 


d, în timp ce în alte situaţii ab 


dacă ae i 0) 
500 


a = bi; „a este aproximativ egal cu d* şi înseamnă că |a—d| 4 cr 
unde c este ales de contextul problemei. 


Y $ 2. Surse de erori 


„Orice rezultat al unui calcul numeric poate fi influențat (uneori dezas- 
truos!) de diferite tipuri de erori. 


) Erori în datele de intrare. De exemplu dacă coeficienții unui sistem 
liniar sînt fracţii ordinare, la transformarea lor în fracții zecimale se introduc 


erori (£ = 0.3333333). ANt exemplu, considerarea în calculul cu numere ira- 


_ ţionale (7, VZ) a unui număr finit de zecimale. 

Uneori datele de intrare sînt rezultatele unor măsurători care sînt siste- 
matic afectate de erori. 
„În problemele tehnico-ştiinţifice (rezolvarea cu calculatorul a ecuaţiilor 
diferenţiale, de exemplu) apar foarte des sisteme liniare algebrice în care 


E: y 5 


coeficienţii sînt integrale ce trebuie calculate aproximativ şi deci coeficienții 
sistemului liniar sînt afectați de erori. 

b) Erori de rotunjire în timpul calculului. Dacă facem un calcul utilizînd 
numai £ cifre în toate numerele cu care lucrăm, atunci produsul exact a două 
numere este format cu 24 sau 2/ — 1 cifre. Dar noi.nu putem utiliza la pasul 
următor de calcul decit rotunjirea la £ cifre. Astfel sc produc erori care pe pat- 
cursul unui lung proces de calcul pot afecta serios rezultatele. 

c) Erori de trunchiere. Acestea sînt erori ce provin din limitarea la o „parte 
din termenii unui proces de calcul înainte de a ajunge la limită. De exomplu, 
Pi Pe 4 aa 3 1 Y 3-0 
din suma progresiei geometrice infinite 1 + «pe +3 i rețibână; 
numai primii patru termeni, diferența dintre suma progresiei şi acest rezultat 

se numește eroare de trunchiere. 

Tot o eroare de trunchiere facem şi atunci cînd aproximăm derivata unei 
funcţii f(x) într-un punct +9 prin raportul: 


Po + D= 102), 
h 


unde H este „mic“. 
De asemenea acelaşi tip de eroare se face cînd pe o porțiune mică în jurul 
unui punct aproximăm o funcție cu tangenta la graficul său în punctul ri 
pectiv. 
d) Erori de simplificare a modelului şi erori comise de om şi maşină. 
Uneori, în practică, în unele modele, idealizăm anumiţi parametri pentru 
a putea prelucra modelul. 
De asemenea se poate greși în programul pentru calculator sau în perto- 
rarea cartelelor, erori de operare ctc. 


$ 3. Erori relative şi erori absolute 


Fie ă o valoare aproximativă pentru cantitatea a (a este valoarea exactă) 
avem: 

Eroarea absolută în ă este dă —a. 

Eroarea relativă în ă este (3 — a)/a dacă aş 0. 

Foarte des eroarea relativă se dă în procente. De exemplu, 4%, înseamnă 
o eroare relativă de 0,04. 

Acum „a — ă“ este corecţia ce trebuie adăugată lui ă pentru a obțin 
valoarea exactă. 

Corecţia și eroarea au aceeași valoare absolută dar semne contrarii. 

Eroarea poate fi pozitivă sau negativă dar noi putem introduce o mar- 
gine pentru eroare care este întotdeauna un număr pozitiv. 

Dacă a = ă + e înseamnă că |ă —a| se 

De exemplu a — 0,462 + 0,011 înseamnă: 


0,451 < a < 0,473 


$ 4. Sisteme de numerație. Reprezentarea numerelor 
în virgulă mobilă 


_ În general sîntem obișnuiți să gîndim şi să calculăm cu numere în baza 10, 
el numărul 257,424 reprezintă: 


2 - 10224 5 - 1014-77 - 1004 4-10-1 4+2-10-2+4+4-10-8 
Orice număr real în baza 10 are o reprezentare unică în forma de mai 


Se exceptează cazul în care, în scrierea pozițională a unui număr zecimal” 
vaza 10) cifra 9 apare de o infinitate de ori. 


Exemple: 

1) numărul 99,9999 ... reprezintă acelaşi număr ca 100 

2) numărul 2,219999 ... reprezintă acelaşi mumăr ca 2,22. 

În general pentru reprezentarea unui număr real putem utiliza o bază 


recare B (număr naturaj) cu B>2 şi atunci se poate arăta (cu excepții ca 
sus) că orice număr real pozitiv admite o reprezentare unică de forma: 


ap B+ apa BI uita B+ao rau" Bitrăz:B* +... 


nde coeficienţii „a“ se numesc cifre în sistemul cu baza B şi sînt numere 
egi pozitive cu proprietatea că: 


> 0osa,<B-—1. 


vantajul mare al scrierii unui număr într-o bază care conduce la o scriere 
oziţională de forma: 


Amtm-a see los da da e: 


unde ay-a +-: aaa reprezintă partea întreagă iar O, a_. a_a ... Partea /racționară 
„a numărului real, este că putem da reguli simple pentru operaţiile aritmetice. 
Cu cît baza este un număr natural mai mic cu atit aceste reguli devin 
mai simple. Acest argument stă la baza faptului că cele mai multe calcula- 


toare opercază în baza 2 (în sistem binar cu două cifre O şi 1). 
Amintim regulile de adunare şi înmulțire în sistemul binar: 


0+0=0 0-0=0 
o+i=1 0-10 
1+0=1 1-0=0 
1+1=10 ici. 


Alte sisteme de numerație utilizate sînt cele octale cu cifre de la O la 7 
şi Hexadecimale cu cifre de la 0 la 9 şi literele A, B, C, D, E, F de la zece 
la cincisprezece. 

De exemplu: 


(13.25) = (1101,01)2 = (15,2)8 = (Die 


unde prin Oa (e: Os: (is am notat baza de lucru. 


Virgulele utilizate aici separă partea întreagă de partea fracționară şi 
se numesc respectiv: virgulă zecimală, virgulă binară, virgulă octală, virgulă 
hexadecimală. 

Un calculator este echipat de regulă cu două tipuri de operaţii. aritme- 
tice, anume calculul cu virgulă fixă şi calculul cu virgulă mobilă. Virgula 
este zecimală dacă baza este 70, binară dacă baza este 2 ctc. 

In cele ce urmează ne vom ocupa în special de baza zecimală. 

Calculul cu virgula mobilă zecimală se efectuează cu un număr constant 

„de cifre după cum se va observa în cele ce urmează. 

Ne oprim acum mai îndeaproape la calculul cu numere reprezentate în 
virgulă mobilă. 


Numim reprezentare: zecimală normalizată cu virgulă mobilă a numărului 
a o reprezentare de forma: 


a= m: 10, unde 0.1 <im|<i 


m este număr real pozitiv sau negativ, 
q este număr întreg pozitiv sau negativ. 


O asttel de reprezentare este totdeauna posibilă, pentru a 7 0. 
Variabila m poartă numele de mantisă și este o parte fracţionară, iar g 
se numește exponent și este întreg. 
De exemplu numărul 542,36 va fi reprezentat ca: 
0,54236 - 103, deci 0,54326 este mantisa iar 3 este exponentul. 
Alte exemple: 
1) 4527 = 0,4527 - 10* mântisa este 0,4527 şi exponentul 4 
> 2) 0,0002617 = 0,2617 - 10-2 mantisa este 0,2617 şi exponentul —3 
3) — 42,15 
4) — 0,012 = — 0,12 - 10-1 mantisa este — 0,12 și exponentul —1 


—0,4215 - 102 mantisa este — 0,4215 şi exponentul 2 


Un număr zecimal în virgulă mobilă este memorat, în multe calculatoare, 
cu mantisa şi exponentul în aceeași celulă de'memorie. 

Celula de memorie însă nu are prevăzută decît o singură poziţie binară 
pentru semn și astfel se ridică problema memorării cînd și mantisa și expo- 
nentul sînt negativi. Soluţia este în a aduna exponentului o constantă pozi- 
tivă de care se ține scamă în calculele aritmetice. 

De exemplu, înțr-o mașină zecimală cu două cifre pentru exponent rom 
aduna 50. la toți exponenţii, (ceea ce înseamnă că exponenții din exemplele de mai 
sus vor fi reprezentaţi în interiorul calculatorului ca 54, 47, 32, 49). 

Pe dou ă poziţii zecimale rezervate exponentului putem memora 100 numere 
diferite (numerele de la O la 99). Cum exponentul poate lua și valori negative 
vom considera un interval de valori de lungime 100 simetric în jurul originii: 


[— 50, +49] 
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Deci exponentul real va lua valori în intervalul [— 50, +49] însă în 
calculator va fi memorat (pe cele două poziţii rezervate lui) un număr cu 50 


mai mare. 
într-un calculator numărul de cifre pentru g și pe este limitat. 


Să studiem acum operaţiile în virgulă mobi 


1) Să presupunem că dorim să adunăm numerele 
va = 112,5 în virgulă mobilă memorat ca 0.1125 - 102 


Xa = 15 în virgulă mobilă memorat ca 0,15 10? 

Aşa cum sînt reprezentate numerele nu putem aduna mantisele deoarece 
virgula numerelor nu este aliniată. 

Diferenţa dintre cei doi exponenți este 1, vom deplasa atunci mantisa 
cclui de al doilea număr (numărul cu exponentul mai mic) spre dreapta cu 
un loc, adică scriem: 

Xa = 04015 - 100 
și acum putem aduna: 
3 a = (0,1125++ 0,015) - 10% = 0,1275 - 109 

Observăm că am lucrat aici cu o mantisă formată cu 4 cifre. 

2) Să considerăm un alt exemplu: - 


x Xa = 0,1213 + 1014 0,4691 - 1071. 
Deplasăim din nou mantisa numărului cu expofent mai mic spre dreapta 
cu un număr de poziţii egal cu diferența exponenților, deci în acest caz cu 2 
şi obținem: 
21 a = 0,1213 - 104 0,004691 - 101 = 0,1259 - 101. 
Se observă că am pierdut cifrele 9 şi 1 ale mantisei deplasate în afara 


capacității de calcul (deoarece am considerat mantisa ca avînd patru cifre). 
3) Presupunăm acum că avem de efectuat următoarea adunare: 


0,9917 - 105+4-0,1230 - 105. 


Numerele de adunat au același ordin de mărime (acelaşi exponent) deci 
singura operaţie ce urmează a se efectua este adunarea mantiselor: 


0.9917 + 0,1230 = 1,1147. 


Rezultatul adunării mantiselor este un număr mai mare ca 1. 
“Noi am definit însă mantisa ca o fracţie zecimală subunitară, 
01 <im<1. 


Operația care urmează a se efectua este deplasarea mantisei spre dreapta 
cu o poziţie zecimală și adunarea la exponent a unei unități. Această operație 
se numeşte renormalizare. 


Deci: 
0,9917 : 1034 0,1230 - 102 — 0,1114- 10% 
(din nou prin deplasare spre dreapta am pierdut cifra 7 deoarece calculul 
se efectuează cu 4 zecimale). ' 
4) Fie acum de efectuat produsul a două mumere x, Şi xa în virgulă mobilă: 
0,1763 - 10-2 x 0,2036 - 10%, 


Normal, se înmulțesc mantisele, 
ma * Ma = 0,1763 x 0,2036 = 0,0358, 


iar exponenţii se adună gi + ge = 3. 
Deci, 
Xe Xe = 0,0358 - 10%, 


Nici de data aceasta mantisa (în modul) nu este cuprinsă în intervalul 
[0,1; 1). 

Operația de renormalizare constă în acest caz în deplasarea spre stinga 
a mantisei 0,0358 cu o poziţie zecimală și scăderea unei unităţi din exponent 
(adică x, - x = 0,3580 - 102). 

Observăm deci că operaţia de renormalizare, în urma efectuării unei ope- 
raţii aritmetice consta în „aducerea: mantisei în intervalul [0/1; 1), cu „ajăs 
tarea“ corespunzătoare a exponentului. 


Observaţie. În majoritatea calculatoarelor înmulțirea se efee- 
tuează pe regiştru de memorie dublu (adică o celulă dublă de 
memorie). Dacă se înmulțesc două numere de lungime „t“, rezul- 
tatul de „22“ sau „2/ — 1* cifre este obţinut în acest registru de 
lungime dublă. Dacă după efectuarea înmulțirii este necesa 
renormalizarea (numai cu deplasarea spre stînga în acest caz) 
cifra de pe ultima poziţie (poziţia „") a numărului va proveni din 
a doua parte a registrului dublu de memorie şi va fi deci corectă. 
În exemplul de mai sus au - a — 0,03588468 - 102 în regi 
de lungime dublă iar rezultatul renormalizat pe patru zecimale 
va fi 0,3588 - 102. 


Dăm în figura I.1 organigrama adunării a două numere în virgulă mo- 
bilă. 

Amintim că, în organigrama din figura 1, un dreptunghi reprezintă un bloc 
de calcul și de atribuire (se efectuează operaţiile din membrul drept al sem- 
nului egal şi valoarea obţinută se atribuie parametrului din partea stîngă 
a semnului egal) iar un romb reprezintă un bloc de comparație şi de decizie. 

Presupunem acum că avem un calculator în care sînt prevăzute 5 poziţii 
zecimale pentru mantisă și două poziții zecimale pentru exponent. 

Cel mai mare număr normalizat ce poate fi memorat cu aceste facilități 
va fi: 

0,99999 - 1022 în realitate 099999 - 101 deoarece am convenit să adunăm 
exponentului real numărul 50. 

Iar cel mai mic număr normalizat (în modul) va fi: 


0,10000 - 10% în realitate 0,10000 - 10-50 
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Fig. 1. 1. Organigrama adunării a două numere în virgulă mobilă. 


(an m Bi, a = ma Bt, au xa = ma Br). 


Cînd printr-o operaţie aritmetică între doi operanzi am depăşit una sau 
Ri irniteie, atat: sua sus, (un număr mai mare de 1,0 - 104 — (depă- 
flotantă superioară) sau un număr mai mic de 1,0 - 10-% — (depășire 

antă inferioară) calculul se oprește. 


10,12641 - 10-42 x 0,98763.. 102 
0,45631 - 10%2 x 0,18764 - 1015 


.. Rotunjirea şi rotunjirea prin tăiere (trunchiere) 


a două căi de rotunjire a numerelor la un număr de zecimale („*) 
anume, rotunjirea prin eliminarea tuturor zecimalelor care se găsesc 
apta. poziției £ (rotunjire prin tăiere) şi a doua cale este rotunjirea obiş- 

că printre numerele ce se pot forma cu „1“ zecimale alegem numărul 
i apropiat de numărul dat. . 


Exemple 


Rotunjirea şi pohunjirea prin tăiere (trunchiere ) la trei zecimale în muemerele 
de mai jos se face după cum urmează: 


- 1) 0,1597 se rotunjeşte Za 0,160; se votunjeşte prin tăiere la 0,159 
2) 0.2456 se rotunjeşte la 0,246; se rotunjeşte prin tăiere la 0,245 
3) 0.1574 se rotunjeşte la 0.157; se rotunjeşte prin tăiere la O, 


4) 0,23652 se rotunjeşte la 0,237; se rotunjeşte prin tăiere la 0,236 

Un mare număr de calculatoare utilizează rotunjirea prin tăiere în rezul- 
tatele obţinute după fiecare operaţie aritmetică, însă dat fiind faptul că numă- 
rul de cifre utilizate în calcul este mult mai mare decît numărul de cifre în da- 
tele problemei nu se ivesc neplăceri. 

Aflarea marginii erorii relative pentru un număr aproximat prin tăiere. 
Să presupunem că avem la dispoziţie un calculator cu 4 cifre disponibile pen- 
tru mantisă. 

Să adunăm x, = 0,1142: 10! cu xs — 0,2190. 101 

Obţinem: 

13 = au Xa = 0,1142 - 101 0,002190 - 102 = 0,11639 - 101 

Acest număr se poate exprima sub forma: 

= 0,1163 - 101 40,9: 10-23, 

Exponentul celui de al doilea termen este mai mic cu 4 decît exponen- 
tul primului termen. Întotdeauna cînd lucrăm cu o mașină cu patru cifre în 
mantisă, exponentul celui de al doilea termen trebuie să fie mai mic cu patru 
decît exponentul primului termen. 

Deci dacă avem un calculator cu „£* cifre în mantisa unui număr repre- 
zentat în virgulă mobilă atunci orice rezultat al unui calcul (obţinut prin 
aplicarea oricărei din cele patru operaţii aritmetice) se poate reprezenta (evi- 
„dent înainte de rotunjire) ca: 

B = my * LO m = 10% 
unde m, este un număr format cu „t* cifre. 

Valorile posibile pe care le poate lua m, sînt: 

01 <Iml<i 


în timp ce: 
0s<|ml<1. 

Spunem că facem o aproximaţie prin tăiere (trunchiere) dacă, pur şi 

simplu nu-l luăm în considerare pe mi, şi acesta este cazul practic cel mai des 


întîlnit la calculatoare. 
Eroarea velativă maximă se obține cînd my este cel mai mare și my cel mai 


mic. >, 
Valoarea maximă posibilă a lui mp este mai mică decit 1 iar valoarea mi- 


nimă a lui m, este O,l. 
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JA 3 La "a. 2 AN AP cite, 


Du relative va fi: 
1 --10e-+ 
0,1 - 10 
"rotunjire, prin același procedeu, obţinem o eroare relativă dată de 

0,5 - 10e-* : 


a = 0,5101, 
0.1 - 102 


= 103% 


"În general pentru o bază arbitrară B avem adevărat următorul criteriu: 
Presupunem că lucrăm în baza B cu „/* cifre în mantisă (cifra binară 
ne dă semnul numărului nu este socotită). Atunci orice număr real în 


, mașinii definită de: ] 


AI B+ dacă utilizăm rotunjirea 
că 


= + dacă utilizăm tăierea. 


În cele mai multe calculatoare „a“ se află între 10-% şi 10-15, 

„Să presupunem acum că operanzii cu care lucrăm ar fi exact reprezen- 
ți în maşină, Nu este sigur că rezultatul va fi exact reprezentat (să ne gîn- 
doar la faptul că rezultatul înmulțirii a două numere de lungime „/* are 
cifre). 

Vom nota cu /H(x +). fl —y). fh(a- x). /l(a]y) rezultatele opera- 
lor de adunare, scădere, înmulţire și împărţire în virgulă mobilă pe care 
iaşina le memorează în urma efectuării acestor operaţii și aproximări prin 
e sau rotunjire. 

Conform criteriului precedent asupra erorii relative avem 


NOW 0 
(O) 
e 12) este una din operaţiile +, —,.,/ 
Să răspundem acum la întrebarea: a) Care este unitatea mașinii pentru 
i calcul cu virgulă mobilă cu 5 cifre zecimale pentru mantisă ? 
m — [0:35 1074 dacă se utilizează rotunjirea 
1- 10-4 dacă se utilizează tăierea. 


n) Fie acum x, — 0,12347 : 10% un număr rezultat din calcul. 
Care este marginea pentru eroarea absolută în acest calcul dacă s-a utilizat 


oarea absolută e, este mărginită de e, = 1 - 102-0,5- 105 — 0,5: 10-2, 
c) Pentru numerele vezultate din calcul va = 0,12347 - 1015 şi as = 
„20006 - 10:0 erorile absolute vor fi mărginite de: es — 0,5 - 1013 şi eg — 
10% 


i d : 3 cd zi) 


d) Care este marginea pentru eroarea absolută în calculul celor trei numere 
de mai sus dacă s-a Jolosit rotunjirea prin tăiere? 


e. = 0,L+ 101 

a = 0,1» 1024 

€a = 0,1 - 10% 

Putem cere variabilei x, o precizie absolută de 0,1 - 10-57 Dar o precizie 
de 0,1 - 10-47 Evident că nu, deoarece valorile de 0,1 - 10-% și 0,1 - 10-40 


sint mult prea mici în comparaţie cu eroarea cu care calculatorul ne „furni- 
zează“ acest număr (0,5 - 10-2). 


Observaţie. La efectuarea unui calcul cu virgulă mobilă trebuie, 
să avem permanent în vedere că ultima cifră zecimală a mantisei” 
este afectată de erori de rotunjire sau trunchiere. Deci dacă am 
obţinut dintr-un calcul cu virgula mobilă, cu 5 cifre zecimale 
pentru mantisă, numărul x — 0,23483 - 107 marginea pentru croa- 
rea absolută este de ordinul 0,000005 - 107 „5 - 102 în cazul 
cînd se aplică rotunjirea sau 0,00001 - 107 Oi - 10% cînd se 
aplică truchierea şi deci nu putem cere o precizie mai bună decit 
0,5 x 102 respectiv 0,1 - 105. 


Acum vom arăta că asociativitatea în adunarea în virgulă mobilă nu se 
păstrează. 
Presupunem că lucrăm cu o maşină cu șapte cifre în mantisă. Fie: 


a = 0,1234567 - 100 
d = 04711325 - 10% 
c=—d 
avem: 
JUB+ c) = 0, fa +/H6+ c)) = fa) = 0,1234567 - 100 
pe de o parte, iar pe de alta: 
Țla+ 2) = 0,00001234567 - 1044 0,4711325 - 10! — 0,4711448 - 10% 
ȚUfUa+ d) + c) = (0,4711448 — 0,4711325))10% = 0,0000123 - 104 = 
= 0,1230000 - 109. 


$ 6. Propagarea erorilor 


Dacă a, = 2,32 + 0,03 şi xa — 1,24 + 0,02 atunci care este marginea 
erorii pentru x — 19? 

Cea mai mare valoare posibilă a lui x, este 2,35 iar cea mai mică valoare 
posibilă a lui aa este 1,22, deci cea mai mare valoare posibilă a lui e — x 
este 2,35 — 1,22 — 1,13. 

Similar cea mai mică valoare a diferenţei x, — za este 1,03 
Deci: 


1.03 & x — a < 1,13 y 


14 


ac — 2a = 1,08 + 0,05. 
cazul general dacă si = Fa E ea, aa = Xa tt ce 


Ta — ea — (Ba ke) < aa — ae < Xitea — (Xa — ea), 

Da — Xa — (ea tea) < xi — Xa S Xa — Xa (eat ea), 

au — aa = Xa — Fat (eat ea). 

n inducție ajungem la următoarea: 

ulă. La adunare și scădere, marginile pentru eroarea absolută sînt 


suma marginilor erorilor absolute în operanzi. 
orm definiţiei erorii relative (vezi $ 3) avem: 


AA — 7 (unde r este eroarea relativă) 


= sU+5. 
acă Xa şi Xa au erorile relative ri Și ra avem: 


Ta > (ra) xa( Lt ra) = x ra) (Lt ra). 
ează că eroarea relativă în produsul 34% este dată de relația: 


E = — rr) Ur) — 
Xe 


= nt ra re 


că [ri] &1 şi [ra] < 1 atunci putem neglija termenul ya și obținem 
rea relativă în X,%a ca fiind aproximativ 7 + 7z- 
e exemplu dacă ri = 0,03 şi ra = — 0,02, eroarea relativă în produs este: 


TF Tok rura = 0,03 — 0,02 — 0,0006 = 0,0094 0,01. 


procedăm la fel pentru a găsi eroarea relativă la împărțire. 
vem a 


Xa = all) 


Ba = xa(t+- ra) 


Bi ah Ltrn 
Ta a kre 
i na Xa 
E ZI a Lie în 
PA 1472 1+ ra 


20 aici zale le 


Dacă |n] < 1 şi |r] & 1 obţinem: 
adi ta: 
1+re 


Dacă marginile erorilor relative în za, a. sînt gi şi pe atunci pat pa este 
cea mai bună margine superioară pentru |r+ ra] ca și pentru Jr, — va]. 
n concluzie avem următoarea regulă: 


La înmulțire și la împărțire marginile pentru erorile relative în operanzi 
se adună. 

O situaţie des întîlnită în practică, care conduce la rezultate eronate este 
următoarea: 

Să presupunem că avem de făcut o scădere între doi. operanzi, în care dife- | 
zența este mult mai mică decit fiecare dintre ei. 

Fie y = x, — Xa. Să notăm erorile în xi şi +a cu Ax, și Axe atunci: 


lAyl < |Azal +A] 


Ra — Te 


deci 


Aplicaţie: Dacă a = 0,6724 + 3 104 şi 


Xa = 0,6723 E 10-4 


atunci: 


a — ra = 00001 2 0,0001. 
Această anulare a rezultatului poate fi ocolită printr-o scriere convena- 
bilă a formulelor de calcul sau aplicarea unui alt algoritm. 


Să considerăm următorul exemplu: 
ecuaţia x? — 56x+ 1 = 0 are rădăcinile: 


a = 28 — VI83 = 28 — 27,982 — 0,018 + - 103. 


Xa = 28 + VT83 a 55.982 + - 10-3. 


Se observă că deşi rădăcina pătrată s-a calculat cu cinci cifre semniti 
tive, totuşi în prima rădăcină (+) primim numai două cifre semnificative iar 
în xa eroarea relativă este mai mică decît 10-5. 

Cum produsul rădăcinilor -este x,xa = ], putem utiliza următoarea cale 
pentru calculul lui x: E 


Hi 


= 0,01786288 
55,982 


m = 


cu o eroare relativă mai mică decît 10%. 
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Aşadar a — 0,0178629 + 0,0000002 și astfel cu aceeaşi estimare a rădă- 
cinii pătrate obţinem pentru +, cinci cifre semnificative (în loc de două) uti- 
Jizînd informația dată de produsul rădăcinilor unei ecuaţii de gradul al doilea. 

în general dacă |3| < x atunci vom scrie: 


3 
Vu +3+Va 


De asemenea în multe situații este preferabil să utilizăm pentru calculul 


diferenţei: 
cos (+3) — cos x 


produsul, 


— a sin A sin (++ =) 


Dăm acum cîteva reguli simple în vederea executării unui calcul practic 
„mai precis. 
1) Cînd se adună sau se scad numere (un şir de numere) să se înceapă 
cele mai mici. 
11) Să se evite pe cît posibil cu putință scăderea a două numere pozitive 
aproximativ egale. Deseori, după cum am văzut, putem rescrie formula de calcul 
ntru a evita acest lucru. 
III) O expresie de forma a(b — c) poate fi scrisă ab — ac iar (a — b)/e 
aje — bc. Dacă numerele sint aproximativ egale să se facă scăderea înainte 
e înmulţire evitînd astfel erori suplimentare de rotunjire. 
IV) Pe cît posibil să se minimizeze numărul operaţiilor aritmetice. 


Exerciţii 


7. Scrieţi în baza 10 numerele (1000), (100000), (464). 

2. Fie numărul 1/3 (de memorat într-o celulă de memorie cu 4 poziţii 
imale. Cum apare acest număr în calculator? Cu ce eroare? Care este 
ecţia? Același Icru pentru numerele şi V2. 

3. în cazul reprezentării în calculator a numărului 1/3 (indiferent pe 
e poziţii zecimale (în număr finit), coincide eroarea de rotunjire cu eroarea 
„de trunchiere? 

4. Pentru reprezentarea (cu rotunjire prin tăiere) a unui număr irațional 
calculator, coincide eroarea de rotunjire cu cea de trunchiere? 


5. Pentru un calculator ce lucrează cu 7 poziţii zecimale, care este 
itatea ? 

6. Un calculator are rezervate 2 poziţii zecimale pentru exponent. De 

se adaugă, în reprezentarea în calculator, numărul 50 exponentului real? 

in vor fi reprezentaţi pe aceste două poziţii exponenții —7, 6, —22, —11,20? 
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7. Avem un calculator cu 8 poziții binare pentru exponent. Care este 
cel mai mare număr întreg ce „încape“ în cele 8 poziţii binare? Care este 
jumătatea acestui număr? Care este numărul ce se va adăuga exponentului 
real în reprezentarea în virgulă mobilă în calculator? În ce interval se situcază 
exponentul real? 


8. Fie de efectuat următoarele operaţii: 

a) 2,1350 - 7; 8) 0,2567 -e;c) m-e;d) 0,2135—0,212;€) 1,2137 +-22,2137 
pe un calculator cu două poziții zecimale pentru exponent şi 4 poziții zeci- 
male pentru mantisă. 

— Scrieţi numerele în reprezentarea virgulă mobilă (normalizat) cu 
rotunjire. 

— Efectuaţi operaţiile pe celula simplă de memorie. 

— Renormalizaţi. 

— Efectuaţi calculele pe celulă dublă de memorie. 

— Renormalizaţi. 


9. Determinaţi eroarea relativă maximă pentru x, unde 
aa = 2,3 + 0,2, 
aa = 22,3 + 0,2, 
23 = 203,3 + 0,2. 
Care din aceste trei numere are precizia relativă cea mai bună? 
70. Calculaţi eroarea relativă maximă pentru: : 
pri aa, 
23 
unde a, aa, x3 sînt cei de la punctul precedent. 
11. Să se efectueze, cu precizia cea mai bună, adunarea: 
S= mat aa+ ao 
în virgulă mobilă cu 5 zecimale (exponent 2 poziţii zecimale) pentru 
22 = 0,11283, 
22 — 0,000043888, 
23 = 0,0082798. 
72. Calculaţi: 


1 
Wa = — da. 
i. , +5 


pentru n = 0, 1, 2,.... 8, utilizind formula de recurenţă: 


i 
Wat Syma = —- 
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Cum rezultă formula de recurenţă? 
ttilizați trei zecimale în calcul şi rotunjiţi. 
rătați că ya > ya şi ya <0 ceea ce este evident absurd. Explicaţi 


. Utilizaţi în exemplul precedent formula: 
: a UR LI 
Vai ră 3 


iderind că yo 2 y arătaţi că ya 2 0,019 şi calculaţi ve Ve +»: Yo: i 


Capitolul II 
REZOLVAREA SISTEMELOR DE ECUAŢII ALGEBRICE LINIARE 


$ 1. Punerea problemei. Condiţionarea sistemelor 


Să considerăm sistemul de două ecuaţii liniare algebrice cu două, necu- 
noscute: 
e ioal Li d) 
da% ae 


Coeficienţii necunoscutelor şi termenii liberi îi considerăm numere reale. 
Vom nota cu A matricea coeficienţilor necunoscutelor, deci: 


| 


şi să presupunem că determinantul matricei A (îl notăm cu det A) este diferit 
de zero, deci sistemul (1) admite soluție unică. 

În această situație matricea A admite o inversă pe care o vom nota 
cu A: Y 


Au Gu 


A= 


An da 


Avem: 
Az2 Zi 
detă detal 
44 = 
Zu au 
deta detal| 


iar, 


—Aadz3 + aaa | 
det A 
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A4- AA —I7 unde I este matricea unitate (la fel se verifică 
N i 

efiniţie. — Se numeşte lungimea euclidiană a matricei 4 (o vom nota 

1,). numărul pozitiv 


All = (aa PF lasa: Iaz + lam PA. 
culind ||A-2|j obţinem: 
_ (anal? lasat lanl?-+ ass), 


a 
Il Lee 
a . - 
“ae oua 2 leat lan P-tlazt?-F leat? , B 
NA Il = le 2 


Definiţie. — Numărul pozitiv dat de relaţia (2) se numește numărul de 

- opace euclidiană a sistemului (1) (îl vom nota cond,(4)= 
la NA). 

Observaţie. Există şi alte numere de condiționare pentru sis- 

teme liniare. Despre unele vom vorbi în acest context iar despre 

altele vor afla cei ce vor intra în studiul aprofundat al tehnicii 

de calcul. 


Din punct de vedere teoretic sistemul (1) admite soluție unică dacă 
A 4.0. Din punct de vedere al calculului practic lucrurile nu stau 

oc astfel. 2 

Să considerăm următorul sistem, de ecuaţii, iz care ne propunem să-l vezol- 

considevînd coeficienţii săi daţi cu două cifre deci, 


8x + 9y=17 6) 
9x + 10y = 19 


Cu regula lui Cramer rezultă imediat că unica soluție a sistemului este: 
„y= 

Din punct de vedere geometric, a rezolva sistemul (3) revine la a intersec- 

ta două drepte (fig. 11.1). 

Considerind x — 2,121, y = 0, în (3) obţinem: 


8x+ 9y = 16,968 4 
9x-+-10y = 19,089. E 


_ Prin rotunjire, la două cifre semnificative, membrul drept al sistemului 
corespunde cu membrul drept al sistemului (3). În concluzie, soluția 


trucit coeficienții şi termenii liberi ai sistemului inițial au fost dați cu două 
e. 
Din figura dată se vede că dreptele sînt aproape confundate, iar punc- 
2,121, y = 0 nu se găseşte pe nici una din cele două drepte, dar el se 
ește foarte âproape de amîndouă. - 
? 
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= 2,121, y= 0 trebuie acceptată ca fiind la fel de bună ca cea unică: 


x(2,121;0) 


Fig. 1.1] 


Calculind numărul de condiționare al matricei asociate sistemului (3) 
obţinem: 


condy(4) = SI BEI t 100 . 


După cum observăm acest număr este destul de „mare“ în sensul următor: 
dacă efectuăm calculul în rezolvarea sistemului (4) cu trei zecimale exacte, 
atunci obținem soluția « = 2,121, y = 0, pe cînd rotunjind în (4) la doui 
cifre avem soluția: + — 1, y— i. 

Așadar numărul de condiționare ne dă o oarecare imagine despre sensi- 
bilitatea soluției unui sistem liniar algebric la perturbații mici în membrul 
drept sau eventual în coeficienţi. 

În fapt, problema rezolvării sistemului (4) poate fi privită şi în felul, 
următor: la o perturbare foarte mică a membrului drept, vectorul 


326. 


- ŢI 
Ed ee E e | 16.968. 
INCI] [19,098 || 
în soluție avem o perturbare destul de mare căci vectorul 
| |a devenit | eIz1 
E o | 


Să introducem acum pentru sistemul (1) şi un alt număr de condiționare. 
Anume asociem matricei A, numărul pozitiv: 


All = max (la + lan]; Laz + | acel) 


(cu alte cuvinte. valoarea maximă dintre numerele pozitive obţinute prin 
adunarea pe linii a modulelor elementelor matricei A). 


Oăservaţie. — Se remarcă faptul că și acest număr pozitiv 
defineşte o lungime pentru matricea A. 
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De asemenea asociind şi matricei inverse, după aceeași regulă, numărul 
pozitiv ||A-1lla, numărul de condiționare a matricei A (pe care-l vom nota 
cu cond.(4)) este prin definiţie produsul numerelor pozitive Alla şi ||4-4l 
deci: 

All * Alla = conda (4)- 

De exemplu pentru sistemul (3) avem: 

|l4lla = max (7; 19) = 19. 
| 


Cum 4 = | 199|] avem: 
9 —s8 


JA“ = max (1—10]4+]9|; 191 +1—81) = 19. 

Rezultă că: 

cond (44) = 19 - 19 = 361. 

Observăm că și acest număr de condiționare este de același ordin de mă- 
rime ca şi condy (4). 

Metodele de'rezolvare a sistemelor liniare algebrice pe care le dezvoltăm 
aici se justifică prin faptul că metoda lui Cramer de exemplu este foarte 
neeconomicoasă (calculul unui determinant necesită foarte multe operaţii) 
şi practic utilizabilă numai pentru sisteme de două sau trei ecuaţii cu tot atitea 


necunoscute. 
Problemele practice dau naştere la sisteme de ecuaţii de ordinul sutelor 


sau chiar miilor de necunoscute şi atunci au trebuit căutate metode eficiente. 

Vom descrie în cele ce urmează, așa-zisele metode directe (datorate lui 
Gauss) adică metode prin care soluția se obține într-un număr finit de paşi 
(un număr finit de operaţii). 

Apoi ne vom ocupa, pe scurt, de metodele iterative (indirecte) în care 
soluţia sistemului se obține printr-un proces de trecere la limită. 

Evident în metodele iterative nu putem efectua o infinitate de operaţii, 
deci ne oprim după calculul unui număr finit de aproximante (se face în acest 
caz o eroare de trunchiere). 

în metodele directe, dacă s-ar lucra ideal (deci fără erori de rotunji 
lucru imposibil dealtfel, s-ar obţine soluţia exactă. 

în timp ce în metodele directe erorile de rotunjire devin foarte periculoase, 
conducînd la soluţii fără sens dacă sistemul este rău condiționat, în metodele 
iterative, în ciuda erorii de trunchiere erorile de rotunjire nu se acumulează. 


$ 2. Metoda eliminării a lui Gauss 


Pentru ilustrarea acestei metode directe vom utiliza un sistem de trei 
ecuaţii liniare algebrice cu trei necunoscute, despre care presupunem că admite 
soluţie unică: E 


asa + arata + dusa = bi . 
aa1X%4 + Aaaa + da3X3 = ba (5) 
2373 + daoXa + d33%3 = ba- 
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Cel puțin unul din coeficienţii dam. aa. as este diferit de zero. căci altfel 
sistemul (5) n-ar mai fi cu trei necunoscute. Dacă au este egal cu 0 rearan- 
jăm ecuaţiile astfel încît coeficientul lui +, din prima ecuație să nu fie zero. 

Soluţia sistemului rămîne, evident, neschimbată dacă schimbăm ecuaţiile 


între ele (De ce?). 
aa 


Definim un multiplicator pia — 
&u 


Înmulţim prima ecuaţie din sistem cu ma şi o scădem din ecuaţia a doua 
(atenție acum „prima” şi a „doua“ se referă la ecuaţiile deja rearanjate dacă 
a fost cazul). - 

După efectuarea operaţiei de scădere, ecuația a doua din sistem va fi: 

(az — Mad) + (ase — Padua) a + (asa — 2219) Xa = 


= ba — mah. (6) 


Dar: 


a. 
azi — Madu = da — 2 au = 0 
&u 


observă de ce a fost 


deci z, a tost eliminat din ecuaţia a doua (acum s 


a 
ales ma =] 
du, 


Notăm: 


aa — Modus 


Giza — Madas 


da — Maha, 
Şi astfel (6) devine: 
atata + atata = Bi. XE 


Deci a doua ecuaţie din sistemul (5) va fi dată de relaţia (7) 
Soluţia acestui nou sistem este aceeași cu soluţia sistemului iniţial. 


Similar, definim ms = =! și înmulțim prima ecuaţie cu ms şi o scădem 


di A treia, Fluaimita adn ou eoefieieatal Ai, 2 ai bb eta 
ascet adsna = 6) 
unde: 
dis = asa — sasa. 
. aia = dos Palas, 
06 == boy Pisa 
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iza 


ind în sistemul iniţial a doua ecuaţie prin (7) şi a treia ecuaţie 

bținem sistemul: : 

aaa aaaăa + asa x9 = bi 
Aiax tata = bi (9) 
ajaxa ++ aja xa = b. 

temul obținut are acecași soluție ca sistemul inițial și prezintă avan- 

i că a, nu mai apare în ultimele două ecuații. 

trecem acum la eliminarea lui xe din ultimele două ecuaţii. 

ă aie — 0, 'permutăm ultimele două ecuaţii. Dacă și ala şi ala sînt 

u zero, am avea o infinitate de soluţii sau nici una, ceca ce nu este 

il deoarece am „admis că sistemul iniţial are soluție unică. 


gem multiplicatorul m 


acă înmulțim a doua ecuaţie din sistemul (9) cu m; şi o scădem din a 
obținem: 


(32 — mă - ai) * Xa (asa — mg * aia) - aa = by — ms * bi 
nou az — m - ala — 0, şi notînd 

Asa = Aaa — Mg * aa 

= bi — mg 


aia - 29 = 8. (0) 
stfel putem înlocui sistemul (9) prin sistemul echivalent: 
ant t aaa aaa Xa = bi 
| diaXa + aia Xa = bi (1) 
aia = d. 
Matricea asociată coeficienţilor necunoscutelor sistemului (11) este su: 


rior triunghiular 
Pau aa da 
lo “aţa aie 
ig7 0. asa 


U= 


Sistemul (11) se rezolvă foarte simplu, prin substituție, începînd cu 
uaţia a treia. 


Amintim că au 7 0, ala ș 0 iar al este de asemenea diferit de zero 
deoarece în caz contrar ar rezulta că matricea asociată coeficienţilor siste- 
mului (5) ar fi singulară (cu determinantul zero). 

în ccuaţia a doua şi a treia din (11) renotăm coeficienții cu dez az, ba 
xespectiv ass, Ba. Aceste ecuații se scriu: 

+ azaaa = ba 
asa = Ba. 

Deşi prin această renotare nu mai putem identifica pe caz, aaa, ba respectiv 
aaa, ba cu coeficienții ecuației doi şi trei din (5) renotarea este posibilă deoa- 
rece nu mai folosim niciodată acest sistem. Renotarea prezintă avantaj la 
scrierea programului de calcul pentru a nu introduce noi variabile ci folo- 
sirea unor locaţii din program rămase libere. 

Operaţiile folosite în rezolvarea sistemului (5) asupra coeficienţilor şi 
termenului liber nu modifică soluţia sistemului. 

Acum să trecem la descrierea generală a algoritmului lui Gauss. 

Să considerăm sistemul liniar de n ecuaţii cu n necunoscute, de asemenea 
presupunem că sistemul admite soluţie umică, deci că matricea coeficienţilor 
necunoscutelor este nesingulară: 


aa 


Aaa Goa + «ce kr Gay «cs F Ara Xa 
Gai aaa --- + Gajă AF ază, 
ei pe cae E (2) 
Aaaa aaa + «eo + ag» 


Dacă notăm: 


atunci sistemul (12) se poate scrie în forma matriceală: 
. Ax=b. (13) 


Presupunem că a d 0. Eliminăm pe x, din ultimele (n—1) ecuaţii 
scăzînd din ecuaţia „i prima ecuație multiplicată cu: 


a 
ma =, 
du 


deci avem n—1 multiplicatori). j 
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mele a — 1 ecuaţii nu vor mai conţine necunoscuta a; vor deveni: 


e) 


aţa + «af = 22 


a = ay — maay 
39 = dp — mad 


= 2, Bec n. 


„Acest sistem este un sistem de n—1 ecuaţii în necunoscutele aa, Xa, Ape: 


Di 
Dacă aţ) 0 similar putem elimina pe aa din ultimele (n—2) ecuaţii 
e ancscutele Xa, Xa deo ae 

Fie: 


ap; 
OSI. (o PND SICA 
aş 
eficienţii acestui sistem sînt daţi de: 
at) = af — mezat? 
1 = 9 — mal, î 3 3, d, ca ne 
Elementele au. a, af, ... ce apar în timpul procesului de eliminare a 
ecunoscutelor poartă numele de elemente pivot. 
Dacă toate elementele pivot sînt diferite de zero, după n—1 paşi obținem 
singură ecuaţie: 
aţa, = 9. 
Acum, colectînd prima ecuaţie din fiecare pas obţinem: E 


au a af n a x 200 


(14) 
aţa = 

de am notat: a) = a; 4 = d; pentru coeficienţii sistemului iniţial pe 
e dorim să-l rezolvăm prin metoda Gauss. 

Am folosit în (14) o notație mai complicată pentru a arăta exact cum 
transformă, coeficienţii și termenul liber în metoda lui Gauss. Remarcăm 
in nou, ca și în cazul sistemelor de trei ecuații, că la scrierea programului 
e calcul este mai simplu să revenim la vechea scriere a ecuaţiei. 


Sistemul (14) este un sistem triunghiular de forma: 


Gatit Gaza --- + Gaman ina = bi 
Gegăat --- + da-ti, 


FE (5) 
Dacă presupunem că ay 4 0, î= 1, 2, „n atunci necunoscutele x, 
aa-as «+» 4 se determină astfel: 
=: (din ultima ecuaţie) 
“n. (din penultima ecuaţie 
(din p ţie) (6) 
Relaţiile (16) se pot scrie compact astfel: 
d — 35 dai 
> DR a RP (7) 


a, 
S LT) 
Matricea coeficienţilor sistemului (15) se numește superior triunghiulară. 
În cazul în care am avea o matrice inferior triunghiulară, deci cu zerouri 
diasupra diagonalei principale, rezolvarea acestui sistem este acecași numai 
că se vor scoate necunoscutele în ordinea xy, Xa: ---. Xa: 
Dacă studiem formula (17) observăm că pentru rezolvarea unui sistem 
triunghiular se fac „n“ împărțiri şi 
PI 1 [E să, 
Şi D= e mn — D adunări și tomulţiri. 


Dacă „n“ este „mare“ atunci m? > n şi 
1 1 
— -nn— am. 
2 d ) 2 


Evaluăm în calcul numărul de operaţii pentru a ne da seama cît timp 
calculator necesită o anume problemă şi deci cît costă rezolvarea ei. 


$ 3. Formularea matriceală a algoritmului lui Gauss 


Urmărind cu atenţie algoritmul lui Gauss se observă că operaţiile care 
se fac asupra coeficienţilor a, se fac şi asupra componentelor vectorului 
format cu termenii liberi. 
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'ormăm matricea At (extinsa matricii A) din matricea A adăugîndu-i 
oană formată din vectorul B. i 


As=||aa aa se ay dm be (18) 


sînt:, = 
— înmulţirea unei linii cu o constantă, 
_— scăderea unei linii din alta și înlocuirea celei de a doua cu rezultatul 


— rearanjarea liniilor. 
Dacă notăm cu așa pe d atunci matricea A& cu „n“ linii și „n + 1% 
oloane este: 


an ame dim dama 
pna az: aan az mia (9 


mi ma: ma ama 


"În urma operaţiilor descrise în paragraful precedent matricea A! devine: 
Ain Anu 


Aa ami 


3 (20) 
PE PE E ARE Ann 
[i + AR +) a, a, 


mm 


Se observă procedeul de renotare a elementelor matricii A& (vezi orga- 
ama și programul) căci elementele a, din (20) se calculează după cum 
arătat în paragraful precedent. 
în procesul de eliminare se observă că la eliminarea necunoscutei x, ele- 
mentele matricei din coloana £ şi liniile £ + 1, „.-+ m(k — 1, 2, nu a — 1) devin 
nule. 


înmulțim acum ultima linie din matricea (20) cu 


şi obținem: 


(o. 0.0....0;1 Se :) i i A 


Renotînd ultimul element al vectorului (21) cu ap. obținem: 
(0, 0, 0, .... 0, 1. cama): (22) 
Ultima componentă a vectorului (22) este componenta x, a vectorului 


soluţie. 
Înmulțim pe (22) cu a,-a. şi scădem din linia „n—1“ şi astfel obținem: 


(0, 0,0, 


5 En-aun-ae 0 Gp-aamri — Anca nam) (23 


Înmulţind vectorul (23) cu obţinem pe ultima poziţie (deci 


mt 
în coloana n-+l a matricii A&) valoarea necunoscutei x, 


Continuînd procedeul obţinem pe rînd necunoscutele x; 
şi matricea A& devine: 
1 Q «a. 


ami 


au= < 503] (24) 


Oi ră 


nun 


Cpe=8n nea Aa Aa 


Observaţie. — În matricea (24) elementele de pe ultima coloană sînt 
exprimate de relaţiile (17) cu b, înlocuit cu ama. 


$ 4. Organigrama pentru metoda lui Gauss şi programul de calcul 
în FORTRAN 


În paragrafele 2 şi 3 am studiat metoda eliminării a lui Gauss şi s-a arătat 
că pentru fiecare pas „i“ în procesul de eliminare este necesar să avem ay zi 0. 
Dacă ay = 0 se caută în coloana i un element diferit de zero, fie acesta ay 
şi se schimbă linia F cu linia ș. Dacă se folosește renotarea înseamnă că avem 
au 40. 

În cursul acestui capitol ($ 6) se va arăta că pentru a obţine o soluţie 
a sistemului algebric liniar, cît mai apropiată de soluţia exactă, nu e necesar 
numai ca a, + O ci trebuie să avem: 


lau] = max la 
ian 


(elementul pivot să fie cel mai mare, în modul, în coloana î)- 
Algoritmul lui Gauss de rezolvare a sistemelor liniare algebrice descris 
în paragraful 2 poate fi urmărit cu- ușurință în organigrama din figura 11.2. 
Urmează programul de calcul corespunzător, în limbaj FORTRAN în 
care sînt demarcate prin cartele de comentariu blocurile de calcul din 


figura WI. 2. 
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Rezolvă sistemul cu înatricea asociată 
sa coeficienţilor necunoscutelor superior | triunghiu- 
tară 


tipăreşte soluția 


Aflarea elementului maxim | în modul) ai 
1| coloanei „i. Fie. k * linia acestui element, 


Motricea este 
singulară 


Este elementul maxim 
in modul) egal cu zero ? 


= nu . 
în 

Schimbă linia „ka elementului maxim cu 
linia „i” Aceeași operajie cu termenul liber 
(schimbă bi cu bj) 


x 
Imparte elementele liniei „i” şi componenta bi 
WI] cu elementul pivot aji 


Scade linia „i” înmulită cu ay, din.linia 
1V| Scade din bj pe bi-ajj 
joi can 


D] 


Fig. 11.2. Organigrama pentru metoda eliminării a lui Gauss 


PROGRAMUL FORTRAN 


pentru metoda eliminării a lui Gauss 


SUBROUTINE GAUSS (A.BN, 
DIMENSION A(N,N), BN) 


AICI ÎNCEPE PROCESUL DE ELIMINAR! 


aaa 


DO 100 IN 
AMA = AD 
IMAX = 1 


e 
C AFLAREA ELEMENTULUI MANIM DIN COLOANA 1 


DO 101] =1N 
IF(ABS(AMAS) — ABS(A(J.D)) 102, 101, 101 


101 

(zi 

C VERIFICĂ DACĂ DETERMINANTUL ESTE DIFERIT DE ZERO 

(2) î 
IE(ABS(AMAX) — EPS) 103, 104, 104 

-103 WRITE (108,3) 
3 FORMAT(I 

sroe 


SISTEM SINGULAR”) 


PERMUTAREA DE LINII, DACĂ ESTE NECESARĂ 
104 DOR = N 
T = AUNAS, R) 
UMAX, KI 
105. AUR) = TJAMAX 
T = BUMAX) 
D(AMAS BD 
BUD) = TJAMAX 
a 
C URMEAZĂ ELIMINAREA NECUNOSCUTE ŞUI) . 
e 


IE(I—N) 106, 100, 100 
106 OJ] =I1+l 
DO 107 ] = II,N 
XM == A(JD 
DO 108K = LN 


103 ACR) me AIR) AR) SM 

107 B(J) = BC) — BD.SM 

100 CONTINUE 
e, 
C AFLAREA SOLUȚIEI PRIN SUBSTITUȚIE 
[3 

DO 109] 2, N 
KaN=J+1 


KR= RAI 
DO 101 =1,R 


10 B(U) = BI) — A, RR): B(KR) 
109  CONTINU 
RETURN 5 
END 


Exemplu mimeric 
Fie sistemul: 


ma at aa = | 
Xa 1,0001 xa+ 2x3 = 2 
aa 2xa+ 2xa = 1 LI 


pe care dorim să-l rezolvăm prin metoda Gauss. 

La pasul 1 în procesul de eliminare avem: 

— elementul maxim în modul în coloana 1 este 1, deci elementul pivot 
va fi au =, 

” — întrucît elementul pivot este egal cu unitatea, prima ecuație a siste- 
mului este deja normată (împărţită la elementul pivot). 

— eliminarea necunoscutei +, din ecuaţia a doua (se scade prima ecuație 
din a doua. deoarece aa = 1) ne conduce la: 


0,0001 xe +ax3= 1, 


eliminarea necunoscutei x, din ecuaţia a treia ne conduce la: 


xa+ 230. , 
temul de ecuaţii devine: a 
+ aa Xatas= 1 
0 --0,0001 xatxa=t 
[ + Xa a =0. 


La pasul 2 al eliminării avem: 
0,0001; as2 = 1, şi elementul pivot este aa (este element maxim 


Mtz aa = 1 
wat xa = 0 


Ă 0,0001 + za = 1 
— eliminăm necunoscuta x din ecuaţia a treia și avem: 


aut aka 
Xa 23=0 
0,9999 = 1. 


La pasul 3 urmează rezolvarea sistemului cu matrice superior triunghiu- 
(se operează cu trei zecimale) 


E = 

da=0— x = —l 

ami l—(—)—1=1 
"Acelaşi sistem, rezolvat cu calculatorul, a furnizat soluţia: 
în 223 = 14000000; ia 1.000099; 3 = 1,000099 


culul s-a efectuat cu 7 zecimale). 


“$ 5. Strategia pivotării 


uticapii 
_ În descrierea algoritmului lui Gauss noi am presupus (vezi $2) că elemen- 
e af? (elementele pivot) sînt diferite de zero, altfel sistemul este singular. 
Să considerăm următorul sistem de ecuaţii: 


ă 5 Xa Xa+ xa=l 
Xa at 2 = 2 (25) 
Xa 20 23 = 1. 

„Sistemul admite soluţie unică (determinantul matricei coeficienţilor este 


ul) xi = le xp le al 
Lă 


dis 


Aplicînd procedeul lui Gauss, după primul pas primim: 
xs=1 
xa-kxa=0. 

= 0 (coeficientul lui x. în prima ecuaţie (26) — 


(26) 


Deci, observăm că 
vezi $ 2). 

Pentru a putea rezolva acest neajuns permutăm cele două ecuaţii şi 
scriem sistemul (26) ca: 


I] 


Xa as =0 
xa=l 
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şi acum obţinem un sistem deja triunghiular care rezolvat ne conduce la 
soluția căutată. 

Să considerăm acum cazul general. Deci la pasul + am obținut aj —0. 
Atunci măcar unul din elementele aţf, 4/92  ... al) este nenul, altfel matri- 
cea A “ar fi singulară. 

Să notăm acel element cu a 3 0. Vom proceda atunci la schimbarea 
între ele a liniilor & şi r apoi vom continua eliminarea. 

Deci orice sistem nesingular poate fi redus la o formă triunghiulară cu 
ajutorul metodei lui Gauss combinată cu schimbări de linii între ele. 

Să revenim acum la sistemul (25) în care coeficientul lui x. din ecuaţia 
a doua (aa2) îl înlocuim cu 1,0001. Rezolvînd din nou prin eliminare obținem: 


| 1 (23) 
—9 999x3 = —10 000. 


Acum, rezolvînd sistemul (28) ca un sistem triunghiular și utilizînd trei 
zecimale obținem: 


m=0, 1.000. (29) 


Am scris 1,000 pentru a scoate în evidență că lucrăm cu trei zecimale. 
Soluţia exactă cu rotunjire la patru zecimale este: 


Xa, = 1000, xa = 1,0001, axa = 1,0001. 


Acest neajuns provine din faptul că a? a devenit după primul pas al 
eliminării 0,0001 (deci zero cînd lucrăm în trei zecimale), așadar un element 
pivot nu are voic să fie nici apropiat de zero. 


Noi am rezolvat în $ 4 sistemul: 
Xa -t Xa-t- Xa Li 

act 1.000 La +- 233 
Xa 2aa+ 23 


prin schimbarea liniilor 2 și 3 şi am obținut soluţia: 
23 = — aa = Xa = 1,000 


corectă în trei zecimale (am utilizat trei zecimale în calculul soluţiei). 
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plul prezentat ne arată că erorile de rotunjire pot avea efect ca- 
al în metoda Gauss, de aceea foarte des la pasul £ urmăm strategia 
ită strategie de pivotare): 


căutăm cel mai mic număr „7“ pentru care: 
al = max |aţp| &sisn 


nbăm liniile & şi r. 

alte cuvinte, pivotul la fiecare pas este cel mai mare element în modul 

oloana considerată. 

Deci, începem prin a căuta pivotul în prima coloană. Linia întii şi linia 

> conține elementul pivot (cel mai mare element în modul din coloana 
se schimbă între ele. - 

rocesul se repetă cu a doua coloană şi așa mai departe. 

d linia ce conţine pivotul este înmulțită cu un coeficient (înainte de a 

scădea din altă linie) ca este de fapt înmulțită cu un număr subunitar 

zi definiţia multiplicatorilor 7, $ 2) şi astfel crorile de rotunjire în liniile 

ot sînt înmulțite cu un număr mai mic ca 1. 


Observaţie. Procesul eliminării a lui Gauss are loc fără schim- 
bări de linii între ele dacă matricea coeficienţilor necunoscutelor 


A => all 
E: 


are proprietatea că este diagonal dominantă, adică dacă: 


leul > Saul 
deci elementele de pe diagonala principală sînt mai mari (sau ega- 
le) în modul decît suma modulelor elementelor din linia respectivă, 


De exemplu: 
| 2 —i 0 
zi 2 —i 
Îi 20 a pl 


o matrice diagonal dominantă de ordinul trei. 
Exerciţiu. Este matricea de ordinul 7 x 7 diagonal dominantă ? 


(2 —1 o o o o lu) 

—1 2 — o o o o 

o —i 2 —1 o o o 

? o o — 2 —1 o [i] 

- [i] [:] o —i 2 —i o 
i 0 310 0 30 2 2 —1 
2 


o o o 0 o; = 
tfel de matrice se numesc tridiagonale şi apar foarte frecvent în practică. 
numesc matrice rare, avînd multe clemente egale cu zero.) 
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$ 6. Sisteme rău condiţionate şi erorile de rotunjire 


Dacă considerăm sistemul Ax = b unde: 
pseă ja aa a. 
Ga Ga aa || 
și Z soluția sa calculată aproximativ, vom numi vector rezidual 7, vectorul 
dat de: r= b— Az. 


Dacă r= 


aştepta că dacă cele două componente ale vectorului rezidual ar fi mici, 
atunci Z să fie o soluţie suficient de bună pentru sistemul considerat. 
Lucrurile nu stau nici pe departe astfel: 
Tati următorul exemplu foarte sugestiv: 
ie: 


==] 


atunci Z este soluţia exactă, deci în mod natural ne-am 


1,2969 0,8648 | 
0,2161 or4a| 

0,8648 (1 
0,1440 | 


1] 
-] 


Presupunem că: 


| 0,9911 | 
| —0,4870| 
Atunci: 


ri 


0,8642 l-] 1,2969 0,8648|| || 0.9911|| || —10e 
Ta ei ez 0,1441] ||— '4870] 104 


Deci componentele lui r sînt destul de mici, de ordinul 10 otuși 2 cu 
componentele date este foarte departe de soluţia exactă, care este: 


A 
|| 

Sistemul considerat este foarte rău condiționat. După eliminarea lui x, 
obţinem: 
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axa — 00, unde: 


aţi) — 0,1441 — ACASA. 0,8648 — 0.1441 — 0,1440999923 = 107%. 


1,2969 


Deci o mică schimbare în coeficientul 0.1441 provoacă o mare schimbare 
în af), deci în ae. 
Dacă inversăm matricea A, obținem: 
pa ups] 0.1141 —0,8648|| 
| —o.2a61 1,2969 || 
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3), — max (10% - (0,1441 -+ 0,8648), 10% - (0,2161 + 1,2969); = 
Ă = 105 - 1,5130 


Al = 2,1617. 


conda(4) = 2,1617 - 1,5130 - 105 2 3,3 - 108. (32> 


laţia (32) ne arată că sistemul considerat este foarte rău condiționate 
trebui să lucrăm în coeficienţi şi în b cu o precizie mai mare de opt 


Observaţie importantă. Chiar dacă elementele matricei A şi 
a vectorului d sînt date exact, totuşi din reprezentarea în calcu- 
lator. (virgulă mobilă) rezultă o matrice rotunjită A cu elemen- 


tele a, unde: 


- laul 


rezentînd unitatea rotunjită în reprezentarea numerelor în virgulă- 
il pe calculator. În cele mai multe maşini de calcul „u* se află între 


Esi 100. 
7. Metode iterative pentru rezolvarea sistemelor liniare. 
Metoda Gâuss-Seidel 
Să considerăm. următorul sistem de două ecuaţii cu două necunoscule 
3x+ y= 4 
34) 
| x —3y= —2. (30), 


ţia sistemului este x= 1, y= 1. 


Scoatem din prima ecuaţie pe x şi din a dona pe y. Obţinem: 


1 
x — 4 —x) 
Ze) 
(35) 
1 
pe (ez a 
3 za x) 
Vom considera următorul proces iterativ: 
at = 
(36) 
1 
Ma = (20 
E Su ) 


plecînd de la o aproximaţie iniţială (de regulă se pleacă de la x = 0, y 
Deci fie + = 0 prima aproximaţie, atunci din prima ecuaţie a lui 
obţinem: 


ut = - 3 (vezi graficul din figura 11.3). 
Apoi din a doua ecuaţie: 
Ass E (- = == » (figura 11.3). 
Continuind procesul de calcul din (36) scoatem: y 
ud = „2 [= 10 iza, 
3 Şi 27 
w pia 26) so 
& 27, s1 


Observăm că deşi am făcut numai două iterații ne-am apropiat destul 
de mult de soluţia sistemului (se vede de asemenea şi pe figura 3). 
Dacă scriem matricea asociată coeficienţilor necunoscutelor sistemu- 


lui (34): 
3 :] 
a i 

observăm următoarele lucruri: 

a) elementele de pe diagonala principală sînt nenule și mai mari în modul 
decît celelalte elemente. , 

2) lau] => lasal = 3- las] = 3-laal 
deci elementele de pe diagonală domină cu raportul trei celelalte elemente. 

Geometric, observăm că panta primei drepte (corespunzătoare primei 
ecuaţii) este —3 (în modul mai mare ca 1) iar a celei de a doua este mai 
mică decît 1. - 
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cum, să vedem ce se întîmplă dacă permutăm ecuaţiile între ele. Deci: 
E ți e A (37) 
y = 4—3x. 
ornind cu acecași aproximaţie inițială (+ = 0, y = 0), obținem 
Iteraţia E y 
o o o 


(38) 


Fig. IL. 4 


Ce s-a întîmplat? 


Panta primei drepte este mai mică decit 1 iar a celei de a doua este mai 
mare decît 1 (în modul). 

Urmează că pentru stabilirea convergenței va trebui să studiem matricea 
coeficienţilor necunoscutelor sistemului. 

Considerăm pentru aceasta forma generală a unui sistem de două ecuaţii 
cu două necunoscute: 


r (za day > bi - (39) 
Za + aa2y = ba: 


Procesul iterativ este definit prin: 


m = A (6 — ay) 

&u (40) 
mb = i, (Be — ax). 

za 


Notăm: 
A = a — a 


Ay = 9 — MD. 
Din prima ecuaţie (39) şi prirha ecuaţie din (40)' deducem: 
Amt = — 12 Apt-b 
du 
şi similar din celelalte două ecuaţii: 
Am = — 2 Aa), 
3 e 
Din ultimele două relații deducem: 


Am m A Aa A sa-i, (4) 
D ai a 
Analog, 
-a, 
= Me: A 
Gun * ase 
Deci: 
Ayn = [> ) A y-2, 
4u * dan 


Continuînd raționamentul, obținem: 


PAL (sia) Ayt-a = [ea Ji Ayo), (42) 
î aaa) au aa) 
La fel se obţine: . 
ACI =); Azto). (43) 
asa * aaa 


„Relaţiile (42) și (43) ne spun că şirul de iterate (x) converge către x 
(prima componentă a soluţiei sistemului (39) şi (y%) converge către y (a 
loua componentă a soluției sistemului (39)) dacă? 


. Bia aa | a. (44) 
i dna * az 
Condiţia (44) este îndeplinită dacă, de exemplu: 
E za > lasel (45) 
Laz] > laaal 


[lg > laal 0) 
Iaz] > laal- i 
alte cuvinte am ajuns la concluzia că diagonala matricei asociate 
ului domină celelalte elemente. . . 
m descrie în continuare metoda Gauss-Seidel pentru sisteme liniare 
i ecuaţii cu trei necunoscute. 

ie sistemul: 


dau Goa Xa aaa 43 = ba (417) 
data + aaa Xa aaa Xa = ba- = 


resupunem că au £ 0, aaa £ 0; aa 7 0. Urmează că sistemul (17) se 
e serie sub forma: 


[si —F aaa Xa das 3 = bi 


a 1 
us De (Ba — aseita — ara) 
dn 
1 
na = E (Ba — aaa — aza43) (45) 
Aaa 


1 
3 = — (bo — dara — Aaaa). 
a 
Vom nota prima aproximaţie cu 44, 449, 09). 
Calculăm noile componente după cum urmează: 


a = (bn — ass — axa?) 
du 


mp — E (be — azi — ass?) (49) 
zi 


mp0 = 7 (ba — aaa — asr). 
- 5 


Astfel s-a încheiat prima iterație. Acum punctul a, 4), aţ devine 
proximaţie inițială pentru iteraţia a doua (ia locul, în memoria calculato- 
lui, punctului (0, 49, x49)). - 

„În general algoritmul este dat de: 


> mo = AL (n — aa” — asa) 
du 
i = 
at — AL (be — aus? — azaxttb) (30) 
îsi ! 


= A (ea — amr! — aaa). 
3 


Putem trece acum uşor la algoritmul Gauss-Seidel pentru un sistem 
de n ecuaţii cu necunoscute. 
Dacă matricea coeficienţilor este A — |a,l;=îs 


"iar vectorul termenilor liberi este: b = 


ap Sa-i, 


Î=Ta 
= (60) 
pentru î= 1, 2, 3, ..-.n- 


Ca primă aproximaţie se consideră, de obicei, un vector u-dimensional 
cu toate componentele nule. 
Calculul se opreşte cînd: 


max | aţi — 7-0] < e, (61) 


Altfel spus, cînd componentele vectorului soluție la două iterații succe- 
sive devin foarte apropiate înseamnă că ne-am apropiat de soluţie. 

Din discuțiile anterioare se poate deduce că metoda descrisă converge 
cînd: 


laul > laal+ alt + anl (62) 


pentru orice i = 1, 2, 3, .... n şi pentru cel puţin un i avem inegalitate strictă. 
Dezi: 


lau] > laal+ aa + lam. (63) 


De exemplu pentru matricea nxn 


Sie dat O 


= i Vai A aa (64) 
. a ti 
[9) —1 2] 
metoda. Gauss-Seidel converge. 
Exerciţiu. — Explicaţi de cet : | 


| 
Dăm în continuare, organigrama (fig. II. 5) şi programul de calcul pentru 
metoda Gauss-Seidel. Urmează de asemenea un exemplu de calcul. 


42 j 
d 


10 
SUM=SUM-Z7 A(LJ) XI) 
5=a 


SUM =SUM:Z 2 A(1,9)- X019) 
EITi 


lă 


xoKk) 
K= 


x), 
N 


DIMENSION A(N, 


Programul FORTRAN 


pentru metoda Gauss-Seidel 
SUBROUTINE SEIDEL (A, B, 


SUM = SUE -t-A(E, I):3109) 


B(S). XO(N). 3 


e 
& 
E SOLUȚIA SISTEMULUI 
e 
8 ER=0 
I1=1 
EI SUM = 0 
nori 
=I121 
TF(10)2,2,3 
3 DO10 ] = 110 
10 
ET 
12 


SUM 4 A(,I).X0(]) 
(BA) — SUM)/A(I,ID) 
ABS(S (0 — X0(I)) 


EEE — AMAXIER, ERID 
1 + 1 
IE(L—N)4,4,5 


3 TF(ER—EPS)7,7,6 
6 DO 2082 1N 
20 SO(R) = XI) 
co 70 5 
7 RETURN 
END 


Exemplu numeric. 


Să se rezolve, prin metoda Gauss 


A:x=bd, unde: 


4 —i 
4 
0 


o —1 


—1 
—1 


—1 
o 
4 

—1 


= 


Seidel sistemul: 


Dăm mai jos tabelul cu primele 10 iterații în rezolvarea acestui 


VECTORUL XO ESTE APROXIMAȚIA INIȚIALĂ. ÎN XI(N) SE OBȚINE 


sistem: 
Nr. iteraţiei = 23 ai zi 

o o [) o o 

1 0.25 0,3625 0,0625 0,40625 

2 0'40625 0:703123 0,203125 0,476303 

3 01476563 0;738251 0,238281 01494141 

Ei 0494141 0,747070 01247070 0;498533 

3 0,198333 0;749267 0;249267 0;499634 : 
â 0,199634 01749817 0,249817 0;499908 : 
i 0499908 0:749954 0;249054 0,499971 

s 0;499971 01749989 0;249989 0,499094 

9 01499994 0;749997 0,249997 01499998 
10 01199998 0,749999 0.249999 0499999. 


Soluţia exactă este: (0,5; 0.75; 0,25; 0,5) şi în cele 10 iterații s-a atins 
recizia de 1,5 - 10% 


xerciţii 


1. Rezolvarea sistemului de ecuaţii: 


duza dra%a 
: dest aaa = ba 


Lai 


revine la a găsi punctul de intersecţie al celor două drepte date de ecuațiile 
sistemului. Să se arate că unghiul făcut de cele două drepte verifică imegali- 


tatea: 
lets al < - condy (A) (se presupune matricea A nesingulară). 


2. Arătaţi că printre toate matricele 2 x 2, nesingulare, ale căror elemente 
sînt întregi nenegativi și nedepășind numărul 100, matricea 


100 99, 
| se) 
are numărul de condiţionare euclidian maximal. 
3. Să se calculeze cond„(A) pentru matricea asociată sistemului: 
5x — 3,31y = 1,69 
(6 — 3,97y = 2,03. 


Cum se schimbă soluția acestui sistem dacă membrul drept devine: 
253 4 (deci s-a schimbat puţin membrul drept). 


4. Este nevoie la rezolvarea sistemului: Ax = b unde: 


i a —L —l 0, "| 

[Zune 1] 

E A || 27 a 2=luj 
o — N] 


de schimbarea liniilor între ele? Explicaţi răspunsul. 
5. Rezolvaţi. pe calculator, prin metoda lui Gauss sistemul: 


| 9, — 2x24 a = 50 


xi + Sa — 3xa = 18 
— 2, + 20 + 73 = 19 
mtilizaţi în calcul 6 poziţii zecimale. 


6. Scrieţi un program FORTRAN pentru metoda Gauss-Scidel şi rezol- 
vaţi prin accastă metodă sistemul: Ax = b, unde: 


| Asa lol] | 

—i 4 o —i 2 

„=! 0 aj zile 
| o —. —. 4] | 


a) În cîţi paşi obţineţi o bună aproximaţie? (lucraţi cu 5 cifre zecimale). 
5) Comparaţi rezultatele cu rezultatele obţinute prin metoda Gauss. 


7. Să se dea o margine interioară pentru numărul de condiţionare 
cond, (4) unde: 


| 1—i || 
A ilbizii atei ae. Pa e-a0, 
Măi ae ăne) 


8. Evaluaţi perturbarea eventuală a soluțiilor sistemului: 


x—2y=—1 
—2x + 4,01y = 2 


dacă. schimbăm componentele membrului drept cu 0,01. Găsiţi soluția siste- 
mului cu aceeaşi matrice a coeficienţilor și cu membrul drept dat de: 


2:53 | Ela 
201 


9. Arătaţi că schimbările de linii sau de coloane într-o matrice A nu 
schimbă numerele de condiţionare cond, (4); cond (4). 

Exerciţii facultative (cap. II) 

Se dă matricea (de dimensiune „x n) 


2 —i IQisa, 
—! 


a) arătaţi că determinantul lui A este n +1; 
2) rezolvaţi ecuaţia de gradul n: 


2—2A 
— 


0. 
0 + 


utilizînd substituția A = 2 cos 8. 
c) Notind cea mai mare rădăcină a ecuaţiei cu A” şi cea mai mică cu 
Aj) calculaţi: 
A 


lim = 
n-o APD 


d) Arătaţi că inversa matricei A este dată de: 


+1 — pu) 


»<u 
ni 
d = 
REED) e 
n+i 


e) Aplicaţi algoritmul lui Gauss pentru a inversa matricea 4, de dimensiuni 
10x 10, utilizînd în membrul drept matricea unitate. 
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Capitolul Il 
REZOLVAREA ECUAȚIILOR NELINIARE 


$ 1. Introducere. Metode iterative 


Dacă /(a) este o funcţie reală de variabilă reală, atunci soluțiile ecuaţiei 
0 în general nu pot fi exprimate sub o formă convenabilă (o formulă 
emplu) şi în consecinţă, pentru obţinerea lor, se adoptă metode itera- 


Aa) 
de 
tive. 

Pentru a ne face o idee despre ceea ce înseamnă metodă iterativă, în 
acest cadru, vom considera ecuaţia: 


x 


F(x) (D 


unde F(x) este o funcţie (de regulă derivabilă de mai multe ori) ale cărei 
valori 1€ putem calcula pentru x aparținînd unui interval. 
Se pleacă de lu un punct +. aproximaţie iniţială, și calculăm şirul: 


a = Flo); Xa = F(aa)i; axa = Flo), (2) 
Fiecare etapă din calcul: 
Xa = Fa) (3) 


se numește iteraţie. 
Dacă şirul (x, converge către o valoare limită z atunci din continuitatea 

funcţiei F(x) rezultă că: 
lim F(x, 


= Flo) 


deci, trecînd la limită în (3) obţinem: 
a = F(a) [6)) 


adică a satisface ecuația x = F(x). 
O interpretare geometrică a acestui calcul iterativ este dată în figura III.1. 
Se vede că, plecînd din x obţinem F(x0). Apoi ducem o paralelă la axa 
Ox prin punctul (o, F(0)) pînă intersectează prima bisectoare y = x, cobo- 
rîm din acel punct perpendiculara pe OX şi obținem pe za(x, — F(x0)), apoi 
Xa — F(x) şi procedeul continuă. 
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Fig. IUL. 


Dacă urmărim graficul funcției din figura III. observăm că în apropie 
soluţiei pantele tangentelor la graficul funcţiei y = F(x) devin foarte mici 
(pozitive). 

Acest fapt ne conduce la a considera comportarea primei derivate în 
vecinătatea soluţiei, 

Utilizînd teorema de medie a lui Lagrange obţinem: 


unde &, se află între xp-a Și aa: . 

în (5) am utilizat şi relaţia de iterare x = F(x,)- 

Urmează că, dacă |£'(3)| < 1 pentru toți « într-o vecinătate a rădăcinii 
care conţine pe xo și x atunci avem convergența asigurată. 

Acum prezentăm riguros demonstrația următoare: 

Teoremă. Dacă ecuaţia + — F(x) are o rădăcină „a“ și F'(x) există în 
intervalul J = (x; |x — a] < ș) şi satisface condiţia: 

|F'()| < m = 1 atunci pentru orice x e J (deci orice aproximaţie ini- 
ţială în intervalul J) avem: 

4) se J.n=0,1,2, 


3) lim a = a 


c) Singura soluţie a ecuaţiei în intervalul J este z. 

Demonstraţie. (a) Demonstrăm prin inducţie. Fie xo e J, presupunem că - 
2-1 € J, atunci arătăm că şi x, e 

Într-adevăr, conform teoremei lui Lagrange avem: 


— a = Ela) — Flo) = FC) (ma 0 


a 


unde î, e J. 
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De unde: 
lua] SF a) (ama, 9) sm lea al s me j 


deoarece xp e J- 
Urmează că Xe J şi (a) este demonstrată. 
(b) Avem: 


la, — als mlăaa —el 


ontinuînd raționamentul obținem: 


la — a sm F(am-a) — Flo) = me IEC) (ama — 9) sm 
(e J)- 


| De unde, recursiv se obţine: 


li — al sm ixo—al. 


Cum 0O<m=< avem lim m" =0. 


Trecînd la limită în ambii termeni ai inegalităţii precedente obţinem: 


lim ja — a] = 0, deci lim x, = a. 


| (€) Presupunem că ecuaţia + = F(x) ar avea în J şi rădăcina f 7 a- 
Atunci: 


i B= F(B) deci, 
(1) (a— 8); ne J 


a—8= Flo) —F(5) = 
de unde: 
la—B| <m-la—Bl<la—3! 


deci o contradicție şi implicaţia (c) este adevărată. 


x 
$ 2. O metodă rapidă pentru extragerea rădăcinii pătrate 


Ecuația x? = c poate fi scrisă şi sub forma x = F(x) unde: 
P ŞI 


F(x) = (+ =) e>o. (6) 


Soluţia ecuaţiei x — F(x) este z— ve (valoarea limită a algoritmului 
ama = F(,)) (fig. 1.2). 
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PT 


Fig. 1.2 


De data aceasta algoritmul se reduce la șirul inductiv: 


(40) 


calculul rădăcinii pătrate a unui număr. 
De asemenea m tabelul cu iteraţiile în calculul rădăcinii pătrate a lui 
2 plecînd de la x = 1,5. 


Fig. II. 3. Organigrama pentru calculul rădăcinii pătrate. 
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Urmează programul FORTRAN pentru calculul rădăcinii pătrate: 


1 FORMAT (GE 12.5) 
READ (103,1) C.EPS 


+" CISI)/2-0 
Tag XD = EPS) 10, 10, 11 
n SI 
So 2912 
10 WRITE (108,2) C, SII 
2 FORMAT (8%. “RĂDĂCINA PĂTRATĂ A NUMĂRULUI', 
1E 12.5, 'ESTE”, E.L 
sroP 
END 


Dăm în continuare tabelul cu primele trei iterații în calculul rădăcinii 
pătrate din numărul! 2: 


ma s su ifie aa 
ş 3 5 
1 15 1416667 
2 1416667 1414216 
3 1414216 1414214 


(UZ = 1414213562 . 


$ 3. Analiza erorilor 


Calculul valorilor funcției F(x) este afectat de erori (de exemplu de erori 
de rotunjire). 

Vom face, în cele ce urmează, o evaluare a erorii după un număr finit de 
iterații. 


Să notăm cu Zi, Za, Zay cui: 2... şirul valorilor calculate. De asemenea, 
notăm cu 3, eroarea din citeai iu? F(z)- 
Avem: 
Za = F(Z) + 3 n=0, 1,2, n, (8 
Fie a rădăcina ecuaţiei a = F(x), deci: 
a= Fa). (9) 
Scăzînd relaţia (9) din (8) obţinem: 
Zu — a = F(2,) — F(a) + 3 (10) 
Aplicînd teorema lui Lagrange, obţinem: 


FC) (Za — 0) a (1) 


unde £, aparţine intervalului determinat de Z, şi de . 


Pasa — 
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Relaţia (11) se mai scrie: 


U—F"(%)) (Zau — 0) > Fa) (Zn — Za) tf 3 (12) 
Presupunem că |F'(t,)| <m<1 şi [3,| < 8 atunci (12) ne dă: 
1 


3. (13) 


m [i 
[Zsa — ela i a lZma 


1—m 
Primul termen al membrului drept estimează eroarea de trunchiere iar 
al doilea termen estimează eroarea de calcul. 
Eroarea de calcul i depinde de eroarea 3, numai din ultima itera- 
= m 
ţie şi deci z, poate fi considerată o cantitate exactă cînd putem estima pe 3: 
Acest lucru ne conduce la următoarea concluzie: Nu este necesar să cal. 
culfim cu mare precizie primele iterate deoarece rotunjirea în aceste iterații 
nu are efect asupra rezultatului final., 
Se spune că metodele iterative se autocorectează. 
Studiind relaţia (13) se-observă că pentru „* suficient de mare în difez 
A : 3 
renţa |2aș: — a| rolul dominant îl are eroarea de calcul ———_-. Aşa că, dacă 
— m 
după un „n“ suficient de mare vom continua iterațiile, vom observa. că dife 


zenţa |Zaș — a] se comportă ca o înărime de ordinul ă 
m 


Exemplu. Pentru calculul lui Vă plecind de la xa = 1,5 obținem: 


CE = — 14415 (rotunjit la 3 zecimale), 
- 1,4142 (rotunjit la 4 zecimale). 
a = 1.41422 (rotunjit ia 5 zecimale). 
1 2 E, ș 
z= (+ cr 2) — 1,41422 (rotunjit la 5 zecimale). 
3 


Pentru Z aparținînd intervalului generat de Z+ și N2 avem: 


LEGI a = 
(a Ze el FR) e 


2—V2|< : 


2. şi observăm că 


F(zi) > E) pentru Y, cuprins în intervalul [Za V2]- 
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LF'(25)| 1 
la — Za] + - -3. 
IEZI papei 
Etectuînd calculele avem: 
x F'(2) = 0.00009 2 1 - 1053, 
i 


IF) 


|za— V2| < 


= 1,0001, 


Îi 1 Li 
E Pic: Fi da Luni 


3 < 1-10"% (precizia în calcul) 
şi cu aceasta: 
|24— VZ| <0,50005 - 10-19 1,0001 - 10%. 
„Se observă că după 4 iterații precizia depinde numai de eroarea de 
calcul: 


= 1,0001 - 10-%. 


1l—m 


Să presupunem acum că x, este aproximaţia de ordin n a rădăcinii simple, 
a ecuaţiei f(x) = 0. 
Atunci; din teorema lui Lagrange: 
Plan) > (a fue 
unde %, aparţine intervalului generat de =, şi a. 
Dacă |/'(39)| > M. pentru xe J atunci: 


IA). 
i (14) 


“Deci avem o estimare a erorii numai funcţie de valorile lui /(x,) (ceea ce 
este foarte convenabil, de multe ori, dacă cunoaştem comportarea. primei 
derivate). 

Să notăm cu (,) valoarea calculată a funcţiei / în x, şi fie: 


Fra) = Par ara). ăla) < 3: 


Vom putea deci exprima cît de bine este aproximată soluţia în funcţie 
de 3. 
Situaţia cea mai bună ar fi să determinăm un +, astfel ca Tla) =0. în 
acest caz valoarea exactă a funcției /(x) satisface inegalitatea: 


Pa] < 3. 


Presupunînd că derivata întii /'(x) nu variază foarte mult în apropierea 
soluţiei + — a, din (14) deducem: 


lu —al< 


3 
TF! i) 


Urmează că cea mai bună margine a erorii pentru orice metodă este e. 
care poartă numele de precizia ce se poate atinge pentru rădăcina a. 


la — al 2 ca unde «= 
TA 
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Dacă |/'(2)| este mic atunci e, este mare şi în acest caz convergența 
șirului (x, este slabă (spunem că problema este rău condiţionată). 
Dacă privim figura III. 4, 


Fig. IL. 4 


observăm că o mică schimbare în funcţie produce o mare variaţie în evaluarea 
rădăcinii. 

Pentru oprirea practică a procesului iterativ nu se folosește e, (nefiind 
cunoscut) ci următorul criteriu: 


Iara — aa la — ama (5) 
la — mal < 3. 


$ 4. Metoda Newton-Raphson 


Metoda Newton-Raphson este o metodă de rezolvare a: ecuațiilor neli- 
niare bazată pe aproximarea curbei (a) cu tangenta sa în punctul 
(ao: /()), apoi următorul punct x, se obţine intersectînd tangenta la grafic 
în Xo cu axa Ox. Procedeul continuă (fig. III. 5). Scriind ecuaţia tangentei în 
(îo, /(x0)) la graficul funcției y — /(%) obţinem: 

o Pao) = fils) (e — aa), de unde pentru ş —0 avem  (presupunind 


Po) 20 
a Za ar azi (7) 
Deci: 
ze (8) 
Se formează astfel şirul iteratelor x x; dat de algoritmul: 
, e, Fan) 
ia ta Au ha unde h, = — 0: (19) 
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9, 


Fig. III. 5 


Funcţia /(x) se aproximează prin tangenta sa în punctul (x, /(x,)) şi 
asa este abscisa punctului de intersecție a tangentei cu axa Ox. 

Sirul (+) converge (nu dăm demonstraţia aici) către o rădăcină a ecuaţiei 
(a) = 0, dacă această rădăcină este presupusă simplă şi alegerea inițială, 10 
este suficient de apropiată de ca. 

Dăm mai jos un criteriu de convergență, uşor de aplicat. i 

Teoremă. — Presupunem că /'(x) 7.0 şi “ nu-şi schimbă semnul 
în intervalul [a, 5] iar f(a) -/(5) <0. Dacă | e | <b-—a, 

a 
Ș 2 22 a,atunci metoda Newton-Raphson converge pentru orice alegere iniţială 
x0 e [a, d]. 

Pentru înţelegerea acestei teoreme a se vedea figura IIL.6. 

Exerciţiu. Vers/icaţi că teorema precedentă se poate aplica pentru rezolvarea 
ecuației 33 — 1,5 = 0 pe intervalul (1, 2]. . 

Din şirul (x) dat de relaţia (19) noi putem calcula, evident numai un 
număr finit de "termeni — de aceea trebuie să spunem cînd oprim procesul 
de calcul (19) astfel încît să avem o aproximaţie dorită a rădăcinii Căutate. 


Criteriul de oprire (vezi programul de calcul) va fi următorul: 


y 


Fie. III. 6 
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Cînd n, devine mai mic decit eroarea admisă în soluție atunci procesul 
de calcul se oprește şi x, devine soluţie aproximativă a ecuației f(x) = 0. 

Dăm în figura III. 7 organigrama și programul FORTRAN pentru metoda 
Newton-Raphson: 


IN 


Tipăreşte : 
E Ati) 
ien af (xi PTT 


Tipăreşte 
ia Kia a Plita lt) 


Fig. 111.7. Organirama pentru metoda Newton-Raphson 


Programul FORTRAN 
pentru metoda Nezton- Ra phson, 


READ (105, 1) XI, EPS 
1 FORMAT (GE 12.5) 


3 HF = F(SD 
HI = HEJFD(SID) 
Su = SI — HI 


IF(ABS(HD — EPS) 10, 10, 11 
WRITE (108, 2) XI, SII, HE, HI 
GE 12 


STOF 
END 
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FUNCTION F(X) 

= se scrie programul pentru f(=) 
END 
FUNCTION FD(S) 
FD = se scrie programul pentru f'(2) 
RETUR! 
END 


Programul de mai sus a fost utilizat pentru rezolvarea ecuaţiei: 
f(a) = sin x — (2)? = 0. 


S-a pornit cu aproximaţia inițială xp = 1,5, Procesul iterativ a condus 
1a soluție (cu precizia de 1: 10-4) în numai patru paşi. 

Exerciţiu. Rezolvaţi aceeași problemă pe calculator complelind programul 
pentru funcțiile F(X) şi FD(X). Folosiţi în calcul 5 zecimale şi porniţi cu apro- 
ximaţia Xe = 13. 

Tipăriți la Jiecare iterație următoarele valori: 


Nr. iteraţici 20 eri Pi 6.2) h 


xerciţii 
1. Pentru a intui convergența sau divergența unei metode iterative de 
rezolvâre a ecuaţiei x = F(x) faceţi reprezentări grafice pentru cazurile 


—1<Fi(3) <0; 0<F'(a)< 1; Fa) > 1; Fa). 


2. Explicaţi cu cuvintele dumneavoastră de ce putem spune că metoda 
Newton-Raphson este o metodă de „liniarizare” plus iterare. 


3, În câţi paşi se rezolvă ecuaţia + — 3 = 0 prin metoda Newton-Raphson? 
Este importantă alegerea iniţială ? 


4. Arătaţi că ecuaţia: 

a8 — 5a2 — 45 ++ 50 — 0 are o rădăcină cuprinsă între 1 şi 2. 
Calculaţi acea rădăcină prin metoda Newton-Raphson alegînd xp = 1 (calcu- 
lul cu 6 zecimale corecte). 

5. Calculaţi yZ prin metoda lui Newton-Raphson cu 6 zecimale exacte. 

6. Calculaţi Y5 utilizînd metoda Newton-Raphson: 


(calculele cu 6 zecimale corecte). 
7. Găsiţi cea mai mică rădăcină pozitivă a ecuaţiei: 
x—tgx=—0 


(cu șase zecimale). 


8. în calculul rădăcinii a a ecuaţiei x — F(x) cu o eroare mai mică decit 
0,5 - 104 s-a obținut: 


Xa = 0.43789 şi as = 0,43814. 


Se ştie că |F'(x)| < 0,4. Ce număr de iterații sînt necesare pentru a fi 
siguri că am atins acurateţea cerută? : 

9. Scrieţi un program de calcul (subrutină) pentru rezolvarea ecuaţiei 
x = F(&), iterativ, Xa = F(x). cu xo aproximarea iniţială şi oprirea cal- 
culului dată de: . 


unde m este dat de Liza 22. 
|, al 


10. Aplicaţi programul precedent la găsirea rădăcinilor ecuației 
3x — cos =0. 


11.) Rădăcina reală a ecuaţiei 32 = x ++ 4 poate fi scrisă: 
pi = e 
1 i E 
2 + ——v321 2 — ——V32z1 
| sr ără V321 + | Ti: 


utilizaţi această expresie şi calculați cu 4 zecimale corecte pe a. 
5) Calculaţi pe a cu metoda Newton-Raphson tot cu 4 zecimale corecte. 
plecînd de la xo = 2. 
12. Calculaţi cu 4 zecimale soluția nenulă a ecuației 5 


z= 4 er. 


Exerciţii facultative (cap. III) 

1) Fic g(a) o funcţie reală de variabilă reală şi xo e R. Vom spune că 
(x) se comportă bine în xo şi vom scrie g e FB(x0), dacă există două inter- 
Vale In = (ao, x0 + e) şi Ia — (ao — e, xo], e > O astfel încît g(x) este convexă 
sau concavă pe I, şi independent convexă sau concavă pe Îa. 


a) Arătaţi că dacă g'(x0) 7 0 atunci ge FB(x0)- 
5) Aparține funcţia g(x) — a? sin 1/x lui FB(0)? 
c) Fie f:R-R cu derivată continuă pe I = (ap — e, do + e) pentru 
un anumit e arbitrar. Presupunem că: 
() = (aro + D) — /(ato — 1) are proprietatea g e FB(0). 
Fie: 


ay = Dee fab, 


arătaţi că pentru hi > 0 suficient de mic d(1) converge monoton către /'(x0). 
d) Pentru g(?) definit la punctul c) cu derivată de ordinul al doilea con- 
tinuă pe I = (0, e), « > O arbitrar şi g, e! e FB(0) arătaţi c 


3 « sh) 
Fi) = ea — 
; 
este monoton crescătoare pentru 4 suficient de mic. 
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e) în ipotezele de la punctul d) presupunînd în plus că avem un şir de 
numere pozitive monoton descrescător (4) care satisface pentru orice & > 1 
relația: 


arătaţi că: i 
pentru toţi d > N (N un anume număr întreg) avem: 


Idea) — dh) < dh) — d). 
/) Tabclaţi derivatele funcţiilor elementare e7, sin x, tg x, ya pe 10, 1] 


înti-o sută de puncte echidistante utilizînd algoritmul dat de schema logică 
următoare (d(4) este dat la punctul c) schema III. 8. 


da 
hi)=1 20; D04p! 
Mn) =-lag Câh) 
= e 


Fig. IIL.S 


Precizăm că f poate fi luat 4, P este precizia cu care lucraţi, iar dacă 
x =0 porniţi cu ho=001-4-P. 
De asemenea am notat cu: 
DUD = (FXOH) —FXOH)oLOH) 
unde literele mari reprezintă reprezentarea în calcul a literelor mici iar ope- 
rațiile încercuite reprezintă aproximarea maşinii la operațiile aritmetice. 
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Răspunsuri 
„Capitalul 1 

1. 5, 32768. 2660, 

2. a = 0,3333 - 100, eroarea = —- 1/3: 107%, corecția — 1/3 - 1074. 

3. da. 

4. numai cînd în reprezentarea a = m - 10% + my - 10 (vezi $ 5) 

0< Im] <05. 

S. 0,5: 10-* pentru folosirea rotunjirii şi 1 - 10% pentru trunchiere, 

6. 43, 56, 28, 39, 80. 

7. NULINII 10027 0020 1-20 20... + 1-20 23—1— 


= 255; » 235 = 128; deci se adaugă 128 la exponentul real 


ge [— 128, 127]. 
8. a) 0,2135: 10% - 0,3142 - 1051 — 0,0671 - 1002 = 0,6710 - 10%, 


d) 0,2135 - 10% — 0,2120 : 10% = 0,0015 - 1050 == 0,1500 + 104%, $ 
2) 0,1214 - 105 4 0,2221 - 1002 = 0,0121 - 10% 4 0,2221 - 102 — 
-— 0,2342 102. 
9. 0.087: 0,009; 0,0009 a treia. 
10. 0,1. A 


11. uk (ek ra) = 012115, pi 
[ea ok aa) + xa = 0,1214, i 
[sta ke Na) o a = 012114. ţ 


12. a) Avem yu + Sya-a -$ 2 aa = um; S 


4) a = In 6 — In 3 —0,182, 
0,0090, Ă LA 
0,050, p 
0.083, i. r, 
a îi 0 aa Idea) i 


e), Eroarea de rotunjire 5 . 104, se amplifică Ia fiecare pas astfel ca 
obţinem în vs o eroare de * 104 = 0,3125. 


13. pg = 0,019, ys — 0,028, ya = 0,058, 

Ya = 0,021; ya = 0.034, ya — 0,088, 

Va = 0,025, ya = 0.043. yo = 0,182 (corect). 
Deci avem stabilitate numerică. 
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Capitolul II 


1, Se exprimă ctg a în funcţie de parametrii directori ai celor două 
„drepte. 
2. Se foloseşte expresia numărului de condiționare dat de lungimea 
euclidiană a unei matrice. 
3. conda(4) = 10.967. Soluţia sistemului iniţial x — 1; > = 1. Soluţia 
sis stecntui periurbat 3 = —12,9; Y:= —20. Perturbarea soluţiei x 22 
= 139; y— 321. 


4. Nu este nevoie. 
5. Xu = 6,154047, xp — 4,313317, x = 3,240207. 
6. După cinci iterații obţinem: 

x = 0,49854, 

Xa 074927, 

Xa m 0,24927, 

za = 0,49963, 
7. cond, (4) > L5-e. 


Soluţia sistemului inițial x = — 1; y=0. Soluţia sistemului per- 
c SA Tal; 3. i 


Se notează 4 —|| _7 au] şi se calculează ||A-Wj„ care număr 


se înmulțește cu 0,01 şi se deduce că apriori putem avea o modificare de 
0,01 în maximum pe componente. 


Exerciţii facultative. 5) Se obține relaţia de recurenţă 
Au mm 2 c05 0 — Ana — Aaa: 


Capitolul III 


3. Un singur pas. Nu. 
4. f(x) = m — Sat — 45x + 50. Se verifică relația: 


= AD -A2) <0. 
5. Aplicind metoda Newton-Raphson obţinem: x = + * (Za ok za). 
7. 3 — 449341. 
8. Se aplică rezultatele din $ 3. 


Ai SUBROUTINE ITER (FUNC, X0, EPS) 


X2) (520 
SM.OX3—X2)/(1— SM) 
GO To 7 


axa 


WRITE (108,8) X3 
FORMAT (13, E 12.5) 
STOP 

END 


10. Luaţi x 20,317 Și mapa = SS veţi obține rădăcina cuprinsă între 
0,316749 şi 0,316754. 
11. 1,7963. 


12. Considerăm /() 
Metoda Newton ne dă: 


ata > nt ni hp — DOD = (a — 1 ez) — 2-2), 
FE) 
Pentru xp = 0,8 obținem x, — 0,79682 deci hy = — 10-% deci rădăcina 


este a == 0,7968. 
Răspunsuri ia exercițiile facultative (cap. III) 


c) avem: d(h) = SU prin derivare obţinem: 
(a) = gh) h = gh) E) - h — (e) — 00) _ 


zh 2-A 
= o, 04 <h (am aplicat formula lui Lagrange pe 10, 4). 
Bă 


Deoarece g s FB(0), pentru k suficient de mic (d(4))' are exact un singur 
semn apoi din n > ha(/n şi Ha fac parte din acei / pentru care raţionamentul 
de mai sus merge) deducem: h 


a 
d() — d) = 4. d(s)ăs 


de unde rezultă monotonia, apoi 
lim d(H) = /"(a0) pentru / cu derivată continuă. 
20 


d) Deoarece g e FB(0) putem presupune că g este convexă pentru A 
suficient de mic, atunci g este crescătoare pentru 4 suficient de mic, deci: 
Fi) =. 
Prin derivare obținem: 


PU) = erp E pp EW=E0, o asa 


= (3) — "09 — n), task. 


62. BeatriceGloria_personal library 


_FB(0) urmează că g” este monotonă. Cum g este presupus 
. că g” > 0 iar e() Ata ieri flo — 2) ne dă g(0)=0. 
că g”” este monoton crescătoare atunci F'(h) > 0 în timp ce g” este 
oton descrescătoare urmează că F'(h) = 0 căci (0) =0 şi” > 0. 

-) Urmînd punctul c) putem scrie: 


jad — a-l = Și = 


Li 


_ gh) 
Ft) 7 | dn 


Lai) — aul > lea ah) una = 


= 1n (pal leu 


Similar avem: : 
ra) —asdr= Și apajeru) —-A0-| an < 
zi | 


lg'(4) —e(h)/h)| 
2 


< In (fu/hesa) - 


de unde rezultă concluzia cerută. 


Bibliografie 


. Bull Gordon, Computational Methods and Algol, George Harop G. & Co LTD., London, 
“Toronto, Wellington, Sydney, 1966. 

. Dahiquist, G. ŞAke Bjărek: Numerical Methods, Prentince Hall, INC., 1974, 

.Dorn S., W.$ D.D. MeCracken: Metode numerice cu programe în FORTRAN, 
Editura tehnică, Bucureşti, 1976. 

4. Ikramov H.: Recueil de problemes de algebre lintaire, Edition Mir, Moscou, 1977. 


PE 


, 63 


„CUPRINS 


Prefaţă -. 


Capitolul I.. 
ERORI. CONCEPTE DE BAZĂ EN CALCULUL NUMERIC .. 


$ 1. Introducere 
$2. Surse de erori K 
$ 3. Erori relative și erori absolute -. 
şa. „Reprezentarea numerelor — în 


“rirgulă mobilă 
$ 5. Rotunjirea și rotunjirea prin tăiere Curaitetitere) 
$6. Propagarea erorilor - 


Capitolul II -. 


REZOLVAREA SISTEMELOR DE ECUAȚII ALGEBRICE 
LINIARE 


. Formularea matriceală a algoritmului 
Organigrama pentru metoda Ii Gawaa și programat 


steme rău condiționaie şi erorile de rotunjire + + 
$7. Metode iterative pentru rezolvarea sistemelor liniare. 

Metoda Gauss-$ 
Exerciţii. 


Capitolul III .. 
REZOLVAREA ECUAȚIILOR NELINIARE... 


$ 1. Introducere. Metode iterative ..- 
$2. O metodă rapidă pentru extragerea rădăcinii pătrate +... + 
$3. Analiza erorilor - 
$4. Metoda Newton-Raphson .-..--- 
Eerci z 
Răspunsuri +... 
Bibliografie... 


36 


